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Molecular hydrogen (H2) is an important component of galaxies because it 
fuels star formation and accretion onto active galactic nuclei (AGN), the two 

o : 

processes that generate the large infrared luminosities of gas-rich galaxies 1 ' 2 . 

J> ! Observations of spectral-line emission from the tracer molecule carbon monox- 

• 1— 1 . 

ide (CO) are used to probe the properties of this gas. But the lines that have 
been studied in the local Universe — mostly the lower rotational transitions 
of J = 1 — > and J = 2 — > 1 — have hitherto been unobservable in high- 
redshift galaxies. Instead, higher transitions have been used, although the 
densities and temperatures required to excite these higher transitions may not 
be reached by much of the gas. As a result, past observations may have under- 
estimated the total amount of molecular gas by a substantial amount. Here we 
report the discovery of large amounts of low-excitation molecular gas around 
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the infrared-luminous quasar, APM 08279+5255 at z = 3.91, using the two 
lowest excitation lines of 12 CO (J = 1 — > and J = 2 — > 1). The maps confirm 
the presence of hot and dense gas near the nucleus 3 , and reveal an extended 
reservoir of molecular gas with low excitation that is 10 to 100 times more 
massive than the gas traced by higher-excitation observations. This raises the 
possibility that significant amounts of low-excitation molecular gas may lurk 
in the environments of high-redshift (z > 3) galaxies. 

APM 08279+5255 was discovered serendiptously during a search for Galactic halo 
carbon stars 4 . It comprises a strongly lensed broad absorption line quasar at z = 3.9 
with a complex optical spectrum and a bolometric luminosity of ~ 10 14 L (where L is 
the luminosity of the Sun) after correcting for lensing magnification; it displays extremely 
strong rest-frame far-infrared emission 5 suggestive of dust-rich gas undergoing rapid star 
formation. The CO(4— >3) and CO(9^8) lines have been detected 3 suggesting that the 
molecular gas must be warm, extended on a scale of up to 1" in the image plane (only 
160-270 pc in the source plane), with a mass of (1 — 6) x 10 9 M©. 

We have used the National Radio Astronomy Observatory's Very Large Array (VLA) 
to detect and image CO J=1^0 and J r =2->l emission from APM 08279+5255. Emis- 
sion in the CO(l— >0) line (Fig. 1) extends over scales of ~ 7" (~ 30kpc at z = 3.9 
for H = 75 kms" 1 Mpc -1 , q = 0.5), well beyond the gravitationally lensed nucleus 6 
where the emission from the CO(4^3) and CO(9^8) transitions is confined 3 . Maps of 
CO(2— >1) detect the northern part of the extended emission and show it to consist of two 
kinematically and spatially distinct features (Fig. 2). The continuum emission detected at 
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23 GHz emanates from the lensed nucleus and is unresolved (< 0.6") with a flux density of 
0.30 ± 0.05 mJy. At 43 GHz, no continuum emission is detected at a level of a = 0.30 mJy, 
and the emission expected from the hot dust in the nucleus 5 is negligible (< 0.04 mJy). No 
appreciable dust continuum is expected underlying the extra-nuclear CO(2^1) emission 
from the cooler dust/gas in those regions. 

The nuclear CO(l— >0) line emission has a velocity-integrated flux density of 0.150 ± 
0.045 Jy km s -1 and a deconvolved size of ~ 0.6" — 1", confirmed by the CO(2^1) map 
(Fig. 2) which shows this emission extending ~ 0.8" along the direction expected from 
matched-resolution 1.6-fim Hubble Space Telescope (HST) images 6 . This size is similar 
to that of the CO (4-^3) and CO (9— >8) emission 3 , while the velocity-averaged brightness 
temperature ratios, R^ 3 = (4 -»• 3)/(l -> 0) and R 2 ^i = (2 -> 1)/(1 -»• 0), are 1.5 ± 0.5 
and 1.35 ± 0.55, respectively, values perfectly compatible with a hot and dense gas phase 
in the nucleus. 

We have considered in detail, and rejected, gravitational lensing as the cause of the 
extended CO(2— >1) structure. There are no obvious optical counterparts and the CO(2— >1) 
maps demonstrate that the two features to the north and north-east and the quasar nucleus 
do not share a common CO(2— »1) line profile, a necessary characteristic if a single galaxy 
underlies them all. Detailed modeling of the lens 6 ' 7 , based on high-resolution HST and 
Keck data, does not suggest the existence of additional lensed features. We thus conclude 
that the extra-nuclear features are due to companion galaxies near APM 08279+5255. 
Moreover, their large angular distance from any central lens caustics precludes significant 
magnification (employing an elliptical mass distribution to represent the lensing galaxy 6 , 
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which accurately recovers the quasar image configuration, the magnification in this region 
is <1.2) and hence the implied molecular gas masses are large. 

Xco is the so-called CO-to-H 2 conversion factor 8 ' 9 ' 10 , ranging from ~ 1 M (K kms" 1 
pc 2 ) -1 in ultraluminous infrared galaxies 11 (ULIRGs) to ~ 5M (K kms -1 pc 2 ) -1 for 
quiescent spirals. The integrated flux density of the two extra-nuclear features in the 
CO(2^1) maps is 1.15 ± 0.54 Jy kms -1 which (for the aforementioned range of Xqo and 
R2->i ~ 1) yields a mass of M(H 2 ) ~ (0.65 — 3.2) x 10 11 M , this value being merely a 
lower limit if sub-thermal excitation of CO(2^1) (i?2-»i < 1) or sub-solar metallicities 
(Xqo > 5M (K kms -1 pc 2 ) -1 ) are important. 

This molecular gas reservoir is ~ 10—100 times more massive than that deduced for the 
nuclear region 3 and neither CO(4^3), CO(9^8) nor high-resolution 210-GHz continuum 
maps give any hint of it. This mass is comparable to that of an entire spiral galaxy and is 
confined within 2" — 3" or ~ lOkpc (Fig. 2). Similar gas-rich companions with no obvious 
optical counterparts, akin to the recently discovered submillimetre galaxy population 12 , 
have been found in the vicinity of other high-redshift systems 13 ' 14 . They have the tell-tale 
marks of intense (Lfir ~ 1O 13 L ), heavily obscured starbursts in merger configurations, 
converting a large reservoir of molecular gas (~ 10 11 M Q ) into stars. In these systems the 
dynamical mass is comparable to the gas mass, a situation possible for the companion 
galaxies to APM 08279+5255, though only observations with higher spatial resolution and 
wider velocity coverage can address this issue properly. 

The failure of sensitive spectral-line observations to detect any extra-nuclear CO(4^3) 
and CO(9^8) emission implies low values (< 0.15) for the global and Rq^s ratios 
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for the massive reservoirs of molecular gas near APM 08279+5255. This has far-reaching 
consequences: at z > 3, millimetre-wave telescopes can only access CO J + l^J, J > 2 
and hence cannot hope to trace the bulk of the H 2 gas. 

The low excitation of lines such as CO(9^8) is hardly surprising, even for a star- 
burst environment, since the average conditions in the molecular gas will not be sufficient 
for their large-scale excitation. If the physical conditions in Orion A 15 are typical of a 
starburst environment, Rg^s < 0.15 (and in most cases Rq^s < 0.05). Moreover, the 
molecular interstellar medium in intense starburst environments may not be just an en- 
semble of warm and dense Orion-type clouds. Where the most energetic star formation 
takes place, the gas may also have a warm but diffuse phase 11 ' 16 , n(H 2 ) < 10 3 cm -3 , gen- 
erated by the disruptive effects of intense ultraviolet radiation, supernovae and the violent 
dynamics associated with mergers. The CO J + l^J, J > 2 emission from such gas 
will be weak, possibly for the entire starburst if this phase is dominant. Starbursts with 
similar far-infrared and low-J CO luminosities can be dominated by drastically different 
gas phases 17 ' 18 . Such differences are expected to occur during the evolution of a starburst, 
depending on the complex interplay between the accumulation of dense gas that feeds the 
star formation (for example, through a merger) and the cloud disruption caused by the 
products of that process, the massive stars. 

Finally, since CO observations at high redshift usually sample emission averaged over 
several kiloparsecs, the possibility of cold and/or diffuse gas dominating the global CO 
emission should not be discounted. The CO J + l^J, J > 2 transitions may then have 
low excitation and thus be very faint, a state not altered by the enhanced cosmic microwave 
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background 14 . The detection of CO at z > 3 means that significant metal-enrichment has 
taken place; low-excitation gas would then represent a previously metal-enriched part of its 
interstellar medium that is not participating in the observed star-formation episode. Such 
excitation gradients have been found in several nearby infrared-luminous galaxies hosting 
powerful cir cum- nuclear starbursts and AGN 19 . 

APM 08279+5255 can be seen as an extreme case of such an excitation gradient, fur- 
ther exaggerated by gravitational lensing of a particularly hot region of the underlying 
system. Differential lensing seems to be preferentially boosting the far-infrared and CO 
emission of the AGN-heated (rather than the starburst-heated) component of the inter- 
stellar medium. This is strongly suggested by the high temperature of the gas and dust 
(~ 200 K) which is not typical of nearby starbursts 20 or clouds heated by young stars 
(e.g. Orion A, T ~ 40 - 60 K). Also, the ratio L FIR /M(H 2 ) ~ 10 4 L© M" 1 is ~40 times 
higher than the highest values observed in ULIRGs and the Lfir/£i.4GHz ratio is an or- 
der of magnitude larger than in star-forming galaxies. This all suggests a heating source 
other than stars and the AGN is the obvious alternative. Differential magnification of 
AGN-heated gas therefore seems to be responsible for the high luminosity of the CO(9^8) 
line. 

While highly excited regions with CO emission amplified by gravitational lensing 
are good markers of the presence of molecular gas at high redshifts, they may give a 
poor representation of its average physical conditions, particularly of its total mass and 
distribution. At z > 3, this excitation bias becomes more serious since millimetre-wave 
interferometers, the instruments usually used for such searches, can detect only CO J+l — > 
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J, J > 2 whose excitation requirements are n(H 2 ) > 10 4 cm -3 and T > 50 K. These values, 
being typical of the average conditions in star-forming clouds, mark a gradual excitation 
turnover. 

Observations of the CO(l— >0) and CO(2^1) transitions, with their minimal excitation 
requirements, may reveal much larger molecular gas reservoirs at high redshifts. The VLA 
is currently the only instrument capable of sensitive, sub-arcsecond observations of these 
lines and will thus be an important tool for unbiased surveys of metal-enriched H 2 gas 
around objects in the distant Universe 21 . 

Methods. 

Our observations took place during 23-24 April 2000 using the VLA in its C config- 
uration. After all overheads, 20 hours were spent integrating on APM 08279+5255. 19 of 
the 27 antennas were equipped with 43-GHz receivers (T sys = 45-100 K); all were equipped 
with 22-GHz receivers, including eleven new receivers with T sys ~ 40 K (c.f. 100-150 K for 
the remainder). 

We exploited the new fast-switching technique 22 ' 23 , recording data every 3.3 s, with 
30 s on the calibrator and 170 s on the source (210 s in the 22-GHz band); the typical 
slewing overhead to the compact phase calibrator, 0824+558, 3.3° away, was 7 s. Pointing 
accuracy was checked every hour. 

This technique yields diffraction-limited images at the highest VLA operating fre- 
quencies over long baselines. Conventional phase referencing would not have been able to 
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track tropospheric phase variations and self-calibration techniques could not be employed 
because of the low signal-to-noise per baseline. 

We used continuum mode, placing emphasis on sensitivity rather than line profile in- 
formation. A dual-polarisation 50-MHz band was placed as close to the expected line centre 
as could be allowed by correlator limitations: 23.4649 GHz for CO J=1^0 (+91 kms -1 
from the line centre 3 ). Bandpass roll-off limited the effective bandwidth per IF to ~45 MHz 
which, at z = 3.9, corresponds to Av ~ 575 kms -1 at 23 GHz. The remaining IF pair was 
tuned to simultaneously observe the continuum at 23.3649 GHz (+1280 kms -1 ). 

We obtained matching velocity coverage at 43 GHz (CO J=2— >1) by placing two 
contiguous 50-MHz dual-polarisation bands at 46.9399 GHz (—20 kms -1 from the line 
centre). The continuum was observed on a separate occasion with both IF pairs tuned 
to 43.3 GHz. The flux density scale was fixed using 3C286; the uncertainty is ~15% at 
23 GHz and ~20% at 46 GHz. Calibration and reduction of the data was standard in most 
respects and the rms noise in all the maps is similar to the expected theoretical limit. 

We cannot completely rule out the possibility that phase errors — caused by atmo- 
spheric fluctuations too fast to be tracked even by our fast-switching scheme — are the 
cause of the extended CO emission, but we consider it remote. This is borne out by an 
extensive series of tests: inspection of the raw visibilities; separate imaging of left- and 
right-hand polarization maps; imaging of the most phase-stable subset of the data; imag- 
ing of the calibrator source. The agreement between the nuclear CO J = 2 — 1 emission 
and its associated 3.5-cm continuum is also strong testimony to the coherence of the CO 
J = 2 — 1 map; the point-like 23.4-GHz continuum emission, observed simultaneously with 
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the neighbouring line emission, provides similar supporting evidence for the resolved CO 
J = 1 — emission. 
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Figure Captions. 

Figure 1. Molecular gas in and around APM 08279+5255, as traced by CO J=l— >0, and contin- 
uum emission from the active nucleus. Top: Naturally-weighted maps at is bs = 23.4649 GHz, middle: 
Z/ Q bs = 23.3649 GHz, corresponding to CO J=l — >0 (V r est = 115.2712 GHz) and its adjacent continuum, 
respectively, at Z = 3.912. The maps were produced from the Fourier transform of the measured vis- 
ibilities, tapered to enhance low brightness extended emission. The dirty beam, identical in both maps 
with FWHM 1.5" X 1.4", is shown bottom left. Contours are at -3, -2, 2, 3, 4, 5, 6, 8, 10 X a, where 
a = 40 //Jy beam - 1 . Bottom: Naturally-weighted map of CO J=l — >0 with the continuum source 
subtracted. The subtraction was performed in the visibility domain and the resulting UV data were then 
Fourier transformed, CLEANed and convolved to a circular beam, FWHM 2.25", shown bottom left. 
Contours this time are at —3, —2, 2, 3, 4, 5, 6 X a. The cross marks the position of lens 6 component 
'A' at RA 08h 31min 41.64s and Dec. +52° 45' 17.5" (J2000). The 0.5" offset between the peak radio 
intensity and 'A' is most likely due to HST astrometric errors. 

Figure 2. Molecular gas in and around APM 08279+5255, as traced by CO J = 2 — 1, overlaid on 
a greyscale image of 8.45-GHz continuum emission from the active nucleus. Top: Naturally-weighted 
map of CO 3=2 — >\ (V re st = 230.5380 GHz) with both IF bands averaged (contours) and overlaid on 
continuum emission at 8.45 GHz (greyscale) at a common resolution (0.5" X 0.5"). Contours are at 
— 3, — 2,2,3,4,5X0", where a = 0.15 mjy beam -1 and the greyscale range is 25(2(7 )-305 //Jy beam -1 . 
CLEAN was applied. Bottom: Naturally-weighted, tapered maps of CO J=2-^l for the two IF bands 
separately, at a common resolution of 0.90" X 0.75" FWHM. Contours are at —3, 3, 4, 5, 6 X O", where 
a = 0.20 mjy beam -1 . CLEAN was not applied. In all cases, the beam FWHM is shown bottom left 
and the cross marks the position of lens component 'A'. 
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Figure 2 
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